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ABSTRACT 

Due to some refinements in the dynamics, we can follow the overall evolution of a 
realistic jet numerically till its bulk velocity being as small as (3c ~ lCP 3 c. We find no 
obvious break in the optical light curve during the relativistic phase itself. However, 
an obvious break does exist at the transition from the relativistic phase to the non- 
relativistic phase, which typically occurs at time t ~ 10 6 — 10 6 5 s (i.e., 10 — 30 d). 
The break is affected by many parameters, such as the electron energy fraction £ e , 
the magnetic energy fraction £g, the initial half opening angle 6q, an d the medium 
number density n. Increase of any of them to a large enough value will make the 
break disappear. Although the break itself is parameter-dependent, afterglows from 
jetted GRB remnants are uniformly characterized by a quick decay during the non- 
relativistic phase, with power law timing index a > 2.1. This is quite different from 
that of isotropic fireballs, and may be of fundamental importance for determining the 
degree of beaming in 7-ray bursts observationally. 
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1 INTRODUCTION 

The discovery of afterglows from gamma-ray bursts (GRBs) 
has opened up a new era in the field. Till the end of Au- 
gust 1999, X-ray, optical, and radio afterglows have been 
observed from about 16, 11, and 5 GRBs respectively (Costa 
et al. 1997; Bloom et al. 1998; Groot et al. 1998; Kulkarni 
et al. 1998, 1999; Harrison et al. 1999; Stanek et al. 1999; 
Fruchter et al. 1999; Galama et al. 1999a). The so called 
fireball model (Goodman 1986; Paczyhski 1986; Meszaros & 
Rees 1992; Rees & Meszaros 1992, 1994; Katz 1994; Sari, 
Narayan & Piran 1996) is strongly favored, which is found 
successful at explaining the major features of GRB after- 
glows (Meszaros & Rees 1997; Vietri 1997; Tavani 1997: 
Waxman 1997a; Wijers, Rees & Meszaros 1997; Sari 1997a 
Sari, Piran & Narayan 1998; Huang et al. 1998a, b, 1999a, b 
Dai & Lu 1998a, b, c; Meszaros, Rees & Wijers 1998; Wijers 
& Galama 1999). However, we are still far from resolving 
the puzzle of GRBs, because their "inner engines" are well 
hidden from direct afterglow observations. 

Most GRBs localized by BeppoSAX have indicated 
isotropic energies of 10 51 — 10 52 ergs, well within the en- 
ergy output from solar-mass compact stellar objects. How- 
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ever, GRB 971214, 980703, 990123, and 990510 have implied 
isotropic gamma-ray releases of 3.0 x 10 53 ergs (0.17 Mqc 2 , 
Kulkarni et al. 1998), 1.0 x 10 53 ergs (0.06 M c 2 , Bloom et 
al. 1998), 3.4 x 10 54 ergs (1.9 M c 2 , Kulkarni et al. 1999; 
Andersen et al. 1999), and 2.9 x 10 53 ergs (0.16 M Q c 2 , Har- 
rison et al. 1999) respectively. Moreover, if really located at 
a redshift of z > 5 as suggested by Reichart et al. (1998), 
GRB 980329 would imply an isotropic gamma-ray energy 
of 5 x 10 54 ergs (2.79 Mqc 2 ). Such enormous energetics has 
forced some theorists to deduce that GRB radiation must 
be highly collimated in these cases, with half opening an- 
gle 9 < 0.2, so that the intrinsic gamma-ray energy could 
be reduced by a factor of 10 2 — 10 3 , and could still come 
from compact stellar objects (Pugliese, Falcke & Biermann 
1999). Obviously, whether GRBs are beamed or not is of 
fundamental importance to our understanding of their na- 
ture. For theorists, the degree of beaming can place severe 
constraints on GRB models. 

How can we tell a jet from an isotropic fireball? Direct 
clues may come from afterglow light curves. As argued by 
Panaitescu & Meszaros (1999) and Kulkarni et al. (1999), 
when the Lorentz factor of the ejecta drops to 7 < 1/8, 
the edge of the jet becomes visible. Thus the light curve will 
steepen by t~ 3 ' 4 , where t is the observed time. This is called 
the edge effect (Meszaros & Rees 1999). Another effect is 



© 1999 RAS 



2 



called the lateral expansion effect. Rhoads (1997, 1999a, b) 
has shown that the lateral expansion (at sound speed) of a 
relativistic jet (7 > 2) will cause the blastwave to decelerate 
more quickly, leading to a sharp break in the afterglow light 
curve. The breaking point is again determined by 7 ~ 1/8. 
The power law decay indices of afterglows from GRB 980326 
and 980519 are anomalously large, a ~ 2.0 (Groot et al. 
1998; Owens et al. 1998; Halpern et al. 1999), and optical 
light curves of GRB 990123 and 990510 even show obvious 
steepening at t > 1 — 2d (Kulkarni et al. 1999; Harrison et 
al. 1999; Castro- Tirado et al. 1999). Recently GRB 970228 
was also reported to have a large index of a ~ 1.73 (Galama 
et al. 1999b). These phenomena have been widely regarded 
as evidence for relativistic jets (Sari, Piran & Halpern 1999). 

However, numerical studies of some other authors 
(Panaitescu & Meszaros 1998; Moderski, Sikora & Bulik 
1999) have shown that due to the increased swept-up matter 
and the time delay of the large angle emission, the sideway 
expansion of the jet does not lead to an obvious dimming 
of the afterglow. Thus there are two opposite conclusions 
about the jet effect: the analytical solution predicts a sharp 
break, while the numerical calculation shows no such sharp 
breaks. The condition is quite confusing. We need to clarify 
this question urgently. 

In a recent paper (Huang et al. 1999c), we have devel- 
oped a refined model to describe the evolution of jetted GRB 
remnants. Due to some crucial refinements in the dynamics, 
we can follow the overall evolution of a realistic jet till its 
expanding velocity is as small as ~ 10 _3 c. Many new results 
were obtained in that paper, e.g., (i) We found no obvious 
break in the optical light curve during the relativistic phase 
itself, i.e. the time determined by 7 ~ 1/8 is not a breaking 
point. But in some cases, obvious breaks does appear at the 
relativistic-Newtonian transition point, (ii) Generally speak- 
ing, the Newtonian phase of jet evolution is characterized by 
a sharp decay of optical afterglows, with the power law tim- 
ing index a > 1.8 — 2.1 . The most interesting finding may 
be that whether the relativistic-Newtonian break appears or 
not depends on £ c , the parameter characterizing the energy 
equipartition between electrons and protons. This has given 
strong hints on the solution to the confusing problem men- 
tioned just above: whether an obvious break appears or not 
may depend on parameters. 

In this paper, we go further to investigate what impact 
will other parameters have on the optical light curves, based 
on the model developed by Huang et al. (1999c). The organi- 
zation of the paper is as follows. For completeness, the model 
is briefly described in Section 2. In Section 3 we investigate 
various parameter effects and present our detailed numerical 
results, mainly in the form of optical light curves. We find 
that the light curve break is really affected by many other 
parameters. Section 4 is our final conclusion, and Section 5 
is a brief discussion. 



time, so that it can trace the evolution of a realistic GRB 
remnant, which is believed to evolve from the highly radia- 
tive regime to the adiabatic one (Dai, Huang & Lu 1999). 
(ii) It takes the lateral expansion of the jet into account. 
The lateral speed is given by a reasonable expression, (iii) 
It also takes many other effects into account, for example, 
the cooling of electrons, and the equal arrival time surfaces. 
The model is very convenient for numerical studies. Here, 
for completeness, we describe the model briefly. For details 
please see Huang et al. (1999c). 

2.1 Dynamics 

Let R be the radial coordinate in the burster frame; t be the 
observer's time; 70 and M c j be the initial Lorentz factor and 
ejecta mass and 8 the half opening angle of the ejecta. The 
burst energy is Eo = 7oM jC 2 . 

The evolution of radius {R), the swept-up mass (m), the 
half opening angle (8) and the Lorentz factor (7) is described 
by (Huang et al. 1999c): 
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where (3 = \J — I/7, n is the number density of surround- 
ing interstellar medium (ISM), m p is the mass of proton, a 
is the co-moving lateral radius of the ejecta (Rhoads 1999a; 
Moderski, Sikora & Bulik 1999), c s is the co-moving sound 
speed, and e is the radiative efficiency. 
A reasonable expression for c s is 

1 1 



c s = 7(7- !)(7- !): 



(5) 



1+7(7-1) 

where 7 ~ (47 + 1)/ (37) is the adiabatic index. In the ultra- 
relativistic limit (7 S> 1,7 ~ 4/3), equation (^) gives c 2 = 
c 2 /3; and in the non-relativistic limit (7 ~ 1,7 ~ 5/3), we 
simply get c 2 = 5/3 2 c 2 /9. 

As usual, we assume that the magnetic energy density 
in the co-moving frame is a fraction £|. of the total thermal 
energy density (Dai, Huang & Lu 1999), B' 2 /8n = f|e', 
and that the shock accelerated electrons carry a fraction 
£ of the proton energy. The minimum Lorentz factor of 
the random motion of electrons in the co-moving frame is 
7c,min = £0(7 - l)ra P (p - 2)/[m c (p - 1)] + 1, where p is the 
index characterizing the power law energy distribution of 
electrons, and m is the electron mass. Then the radiative 
efficiency of the ejecta is (Dai, Huang & Lu 1999) 



2 MODEL 

We use the model developed by Huang et al. (1999c). This 
model has the following advantages: (i) It is applicable to 
both radiative and adiabatic blastwaves, and appropriate 
for both ultra-relativistic and non-relativistic stages. The 
model even allows the radiative efficiency e to evolve with 



(6) 



where t„ x is the co-moving frame expansion time, t' cx — 
R/{-yc), and t' ayn is the synchrotron cooling time, t' syn — 
67rm e c/(a"TS' 2 7 ei min), with ot the Thompson cross section. 
In this paper, we call the jet, whose radiative efficiency 
evolves according to equation (^), a "realistic" one (Dai, 
Huang & Lu 1999). 
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2.2 Synchrotron Radiation 

In the co-moving frame, synchrotron radiation power at fre- 
quency v 1 from electrons is given by (Rybicki & Lightman 
1979) 



V3e 3 B' 
m e c 2 



dK 



F 



(7) 



where e is electron charge, 7 e ,max = 10 8 (B'/1G) _1//2 is the 
maximum Lorentz factor of electrons, dN' c /d"f c is electron 
distribution function, u' c = 37e eB' /(Anm e c), and 



F(x) 



r+oo 

= x K, 

J X 



5 / 3 (fc)dfc, 



(8) 



with K 5 / 3 (k) being the Bessel function. 

In the absence of radiation loss, the distribution of the 
shock accelerated electrons behind the blastwave is usually 
assumed to be a power law function of electron energy, 
dNe/d'fe oc 7,T P . However, radiation loss may play an im- 
portant role in the process. Sari, Piran & Narayan (1998) 
have derived an equation for the critical electron Lorentz 
factor, 7c, above which synchrotron radiation is significant: 
7c = 6irm c c/ (<7t7 B' 2 t). In our model, the actual electron 
distribution is given according to the following cases (Dai, 
Huang & Lu 1999; Huang et al. 1999c): 



(i) For 7c < 7c 
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(ii) For 7c,min < 7c < 7c 



dNL 



7e 



"(P+I) 



(7c,min < 7c < 7c), 
(7c < 7c < 7c, max); 
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dy c I 7, 
(iii) For 7 C > 7, 

dN 
d-y, 

Let Dl be the luminosity distance, be the angle be- 
tween the velocity of emitting material and the line of sight 
and define fx — cosO. Then the observed flux density at 
frequency v from this emitting point is 

In order to calculate the total observed flux density, we 
should integrate over the equal arrival time surface (Wax- 
man 1997b; Sari 1997b; Panaitescu & Meszaros 1998) deter- 
mined by 



=r- 



■01* 



dR = const, 



(13) 



within the jet boundaries (Moderski, Sikora & Bulik 1999; 
Panaitescu & Meszaros 1999). 



3 NUMERICAL RESULTS 

In this section, we follow the evolution of jetted GRB rem- 
nants numerically to see what effects will those parameters 
(such as £ e , #o, #obs, n,p, etc.) have on the optical light 
curves. For convenience, let us define the following initial val- 
ues or parameters as a set of "standard" parameters: initial 
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Figure 1. The evolution of the Lorentz factor 7. The solid line 
corresponds to a jet with "standard" parameters. Other lines are 
drawn with only one parameter altered or one condition changed 
with respect to the solid line. The dashed line corresponds to a 
jet without lateral expansion, i.e., c s = cm/a; the dotted line is 
for 0o = 0.1; and the dash-dotted line is for n = 10 6 cm~ 3 . For 
the meaning of "standard" , see Section 3 in the main text. 
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Figure 2. The evolution of the half opening angle 6. The solid 
line corresponds to a jet with "standard" parameters. Other lines 
are drawn with only one parameter altered. 



energy per solid angle E0/Q0 = 10 54 ergs/47r, 70 = 300 (i.e., 
initial ejecta mass per solid angle M c - s /Q,q w 0.0019M q /47t), 
n = 1 cm" 3 , £l = 0.01, p = 2.5, D L = 1.0 x 10 6 kpc, 
£ e = 0.1, 60 = 0.2, ^ bs = 0, where the observing angle o h s 
is denned as the angle between the line of sight and the jet 
axis. For simplicity, we first assume that the expansion is 
completely adiabatic all the time (i.e. e = 0, we call it an 
"ideal" jet, distinguishing it from the "realistic" jet defined 
in Section 2.1). 

Figure 1 shows the evolution of the Lorentz factor for 
some exemplary jets. In the "standard" case, the ultra- 
relativistic phase lasts only for ~ 10 5 s, and the non- 



© 1999 RAS, MNRAS 000, 000-000 




Figure 3. The effect of the parameter £ c on the optical light 
curve. The solid line corresponds to a jet with "standard" pa- 
rameters. Other lines are drawn with only £ e altered. 



Figure 5. The effect of the parameter 6$ on the optical light 
curve. The solid line corresponds to a jet with "standard" pa- 
rameters. Other lines are drawn with only 8q altered. 




log t (s) log t (s) 



Figure 4. The effect of the parameter £g on the optical light 
curve. The thick solid line corresponds to a jet with "standard" 
parameters. Other lines are drawn with only £j| altered. 



Figure 6. The effect of the parameter O ^ S on the optical light 
curve. The solid line corresponds to a jet with "standard" param- 
eters. Other lines are drawn with only 8 ohs altered. 



relativistic phase begins at about t ~ 10 6 ' 5 s. In short, the 
ejecta will cease to be highly relativistic at time t ~ 10 5 — 
10 s. This again gives strong support to our previous ar- 
gument that we should be careful in discussing the fireball 
evolution under the simple assumption of ultra-relativistic 
limit (Huang et al. 1998a, b, 1999a, b). In the case of a 
dense ISM (n = 10 6 cm" 3 , the dash-dotted line), the ex- 
pansion will become non-relativistic as early as t ~ 10 4 ' 5 s. 
Figure 2 is the evolution of the jet opening angle 9. During 
the ultra-relativistic phase, 9 increases only slightly. But at 
the Newtonian stage, the increase of 9 is very quick. 

The effect of £ e on the optical (R-band) light curves is 
illustrated in Figure 3, from which we can see clearly that: 
(i) In no case could we observe the theoretically predicted 
light curve steepening (with the break point determined by 



7 ~ 1/9) during the relativistic stage itself, this is consis- 
tent with the result of Moderski, Sikora & Bulik (1999). 
Note that in this figure the relativistic phase is restricted 
by t < 10 6 s. Huang et al. (1999c) have given a reasonable 
explanation to this phenomena: the edge and the lateral ex- 
pansion effects begin to take effect when 7 ~ 1/6 (in our 
calculation this occurs at t ~ 10 5 ' 5 s, however the blast- 
wave is already in its mildly relativistic phase at that mo- 
ment and it will become non-relativistic soon after that (i.e., 
when t > 10 6 ' 5 s, see Figure 1). So it is not surprising that 
we could not see any obvious breaks during the relativistic 
phase, they just do not have time to emerge. Another possi- 
ble reason is: at these stages, since 7 is no longer much larger 
than 1, we should be careful in using the analytic approxi- 
mations; (ii) When £ c is small, an obvious break does appear 
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Figure 7. The effect of the parameter n on the optical light 
curves. The solid line corresponds to a jet with "standard" pa- 
rameters, other lines are drawn with only n altered. 



in the light curve, but it is clearly due to the transition from 
the relativistic phase to the non-relativistic phase. The sim- 
ulation by Moderski, Sikora & Bulik (1999) does not show- 
such breaks, because their model is not appropriate for non- 
relativistic expansion (Huang et al. 1999c); (iii) When £ c is 
large, the break disappears. This is not difficult to under- 
stand (Huang et al. 1999c). According to the analysis in the 
ultra-relativistic limit, the time that the light curve peaks 
scales as (Wijers & Galama 1999; Bottcher & Dermer 1999; 
Chevalier & Li 1999) 

/ o\ 4 / 3 

UocH^fj ? c 4/3 (Cl) 1/3 . (14) 

/.From Figure 3 we can also see that with the increase off c , 
tm becomes larger and larger, consistent with equation (|l4|). 
In the case of £ c = 1.0, t m is as large as ~ 10 5 s. Note that 
the expansion has already ceased to be ultra-relativistic at 
that moment. Then we can not see the initial power law 
decay (i.e., with a ~ 1.1) in the relativistic phase, and 
so it is not strange that the break due to the relativistic- 
Newtonian transition does not appear (Huang et al. 1999c); 
(iv) In all cases, the light curves during the non-relativistic 
phase are characterized by quick decays, with a > 2.1. This 
is quite different from isotropic fireballs, whose light curves 
steepen only slightly after entering the non-relativistic phase 
(Wijers, Rees & Meszaros 1997; Huang, Dai & Lu 1998a). 
We thus suggest that the most obvious characteristic of jet 
effects is a sharp decay of afterglows at late stages (with 
a > 2). This will be further proved by other figures followed. 

Figure 4 illustrates the effect of £j| on the optical light 
curves. Again no break appears during the relativistic phase 
itself. Interestingly but not surprisingly, we see that has 
an effect similar to £ e : for small values of £|, there are ob- 
vious breaks at the transition from relativistic stage to non- 
relativistic stage (i.e., at t ~ 10 6 ' 5 s); but for large £b values, 
the break disappears, we could only observe a single steep 
line with a > 2.1. The reason is also similar to that of £ e - 
With the increase of t m becomes larger and larger, as 



Figure 8. The effect of the parameter p on the optical light 
curves. The solid line corresponds to a jet with "standard" pa- 
rameters, other lines are drawn with only p altered. 

can be clearly seen from equation (jlj). When t m is large 
enough (i.e., enters the mildly relativistic zone), the initial 
power law decaying segment (with a ~ 1.1) in the rela- 
tivistic phase will be hidden completely. Then we can only 
see the quick decay in the Newtonian phase. The figure also 
shows that with the decrease of £|, the flux density decreases 
substantially. 

Figure 5 illustrates the effect of Go on the light curves. 
For small Go values, the breaks due to the relativistic- 
Newtonian transition are very obvious, consistent with a 
recent analytic treatment by Wei & Lu (1999). For large 
Go values, the breaks are not striking. The dotted line is 
drawn with Go = 1.3, it in fact can be regarded as an 
isotropic fireball. The most notable difference between the 
solid line and the dotted line is their slope disparity in the 
non-relativistic phase. This is an important difference be- 
tween jet and isotropic fireball, and may be useful in deter- 
mining the degree of beaming. 

The effect of G ODB on the light curves is shown in Figure 
6. The thick solid line is drawn with G OOB = 0. The light curve 
with G bs = #o deviates the thick solid line only slightly. 
But for G OOB > Go, the observed flux decreases severely. For 
#obs = 0.3, the observed peak flux density is dimmer than 
that of ^obs = by 5 mag. But fortunately, their late time 
behavior is very similar. So it is still possible for us to resort 
to systematic deep optical surveys, which are expected to 
find many faint afterglows not associated with any know 
GRBs if GRBs are highly collimated. 

Figure 7 illustrates the effect of ISM number density n. 
According to analysis under ultra-relativistic assumption, 
t m is not related to n (t m tx n°), and the peak flux density 
scales as: Sn.,m oc n 1//2 (Wijers & Galama 1999; Bottcher & 
Dermer 1999). Figure 7 shows these patterns qualitatively. 
Figure 7 also shows that n affects the light curve steepening: 
with the increase of n, the time that the ejecta enters the 
Newtonian stage comes earlier, so that the light curve break 
becomes less and less striking. The effect of dense media on 
afterglows from isotropic fireballs has been discussed by Dai 
& Lu (1999a, b). 



© 1999 RAS, MNRAS 000, 000-000 




Figure 9. Optical light curves under different assumptions. All 
the lines are drawn with "standard" parameters. The solid line 
corresponds to an adiabatic jet, while the dash-dotted line cor- 
responds to a jet without lateral expansion, and the dashed line 
corresponds to a jet that evolves from highly radiative regime to 
adiabatic regime. 



Figure 8 illustrates the effect of p on the light curves. 
From equation (|l4|) we know that t m is related to p. Our 
numerical results reflect the relation correctly. Generally 
speaking, although p affects the slope of the light curve 
notably, it has minor effect on the light curve steepening. 
This again is very different from that of an isotropic fire- 
ball, whose timing index a varies from (3p — 3)/4 in the 
relativistic phase to (15p — 21)/10 in the Newtonian phase 
(Wijers, Rees & Meszaros 1997; Dai & Lu 1999a, b), i.e., the 
increment Aa strongly depends on p: Aa — (27 — 15p)/20. 

Light curves in Figure 9 are drawn under different 
assumptions. The solid line corresponds to a "standard", 
"ideal" jet (adiabatic), the dashed line corresponds to a "re- 
alistic" jet (whose radiative efficiency e evolves according 
to equation ([])), and the dash-dotted line corresponds to 
a jet without lateral expansion (i.e., c s = 0). It is gener- 
ally believed that the lateral expansion effect tends to make 
the light curve steepening more obvious. However, Figure 
9 has just shown the opposite tendency: the steepening of 
the dash-dotted line is obviously more striking than that of 
the solid line. In fact the numerical simulation by Moderski, 
Sikora & Bulik (1999) also shows this tendency. From this 
we conclude that it is the jet edge effect that leads to the ad- 
ditional light curve steepening (with respect to an isotropic 
fireball) at the relativistic-Newtonian transition point. The 
lateral expansion effect tends to reduce the steepening. 

Figure 10 shows the evolution of the observed spectra 
for the "standard" jet. At early times, the spectrum can be 
divided into three segments. But at time of t ~ 10 6 s, a 
fourth segment emerges. It may come from the edge effect 
in Newtonian phase, and it is just this segment that leads 
to the light curve steepening. 



Figure 10. The spectral evolution of a jet with "stan- 
dard" parameters. vS v is in units of ergs cm -2 s —1 . The 
curves from upper-right to lower-left correspond to t = 
10 2 , 10 3 , 10 4 , 10 5 , 10 6 , 10 7 , and 5 X 10 7 s respectively. 



4 CONCLUSION 

We have studied the evolution of jetted GRB remnants 
numerically, following the convenient model developed by 
Huang et al. (1999c). In typical cases, the remnant enters the 
mildly relativistic phase at time t ~ 10 s — 10 6 s (i.e., 1 — 
10 d), and it will become non- relativistic at time t > 10 6 ' 5 s 
(37 d). In analytic approximation, it is usually assumed that 
the expansion is highly relativistic (7 3> 1) all the time. Now 
we should note that these approximations are appropriate 
only for early afterglows (Huang et al. 1998a, b, 1999a, b, 
c). Due to lateral expansion, the half opening angle 6 will in- 
crease with time. But the increment of 8 is very small at the 
ultra-relativistic stage. After entering the Newtonian stage, 
the increase of 8 is very quick. 

It has been predicted that due to the edge effect and the 
lateral expansion effect, the light curve of afterglows from 
jetted GRB remnants should show obvious steepening dur- 
ing the relativistic phase, with the break point determined 
by 7 ~ 1/8. However, our numerical results show clearly 
that no steepening could be observed during the relativis- 
tic stage itself (also see Huang et al. 1999c). But in many 
cases, a striking break does appear in the light curve at later 
stages, which is obviously due to the relativistic-Newtonian 
transition. Typically the power law timing index a varies 
from ~ 1.2 to > 2.1. 

The light curve break due to the relativistic-Newtonian 
transition is affected by many parameters, such as £ e , 80 
and also n. Increase of any of them to a large enough value 
will make the break disappear. We have also demonstrated 
that it is the jet edge effect (Kulkarni et al. 1999; Meezaros 
& Rees 1999) that leads to the light curve break. The lat- 
eral expansion effect just tends to reduce the steepening (see 
Figure 9), this may be completely unexpected to most re- 
searchers. 

Although whether the break appears or not depends 
on parameters, afterglows from jetted GRB remnants are 
uniformly characterized by a quick decay during the non- 
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relativistic phase, with a > 2.1. This is quite different from 
isotropic fireballs, whose light curves steepen only slightly 
after entering the non-relativistic phase. It can be regarded 
as a fundamental characteristic of jets and may be used to 
judge the degree of beaming. 

The effect of p on the light curve break is also very in- 
teresting. For an isotropic fireball, if p is as small as 2.1, 
the light curve will steepen only slightly after entering the 
non-relativistic phase. But for a jet, if other parameters are 
properly assumed, we still can observe a steep break. In prac- 
tical observations, if p can be determined to be near 2 from 
spectral information, and if we observed a steep break in the 
light curve, then possibility is large that we were observing 
afterglows from a jet. 



This means that there are at least two kinds of GRB after- 
glows: one corresponds to an isotropic fireball in a uniform 
medium (i.e., case (a) ), the other corresponds to an isotropic 
fireball in a wind environment or a jet. These kinds of infor- 
mation can provide important clues to our understanding of 
GRBs. For example, at least we know that 7-ray emission 
in some GRBs is isotropic, then the energy crisis is really a 
problem: GRB 971214 and 980703 indicate isotropic energies 
of ~ O.17M c 2 and ~ 0.06M o c 2 respectively. 

Systematic deep optical surveys may provide another 
way for determining the degree of beaming. If GRBs are 
highly collimated, these surveys may reveal many faint 
decaying optical sources. They are afterglows from jet- 
ted GRBs whose 7-ray emission deviates the line of sight 
slightly. 



5 DISCUSSION 



ACKNOWLEDGMENTS 



We have shown that the afterglow from jetted GRB remnant 
is characterized by a steep light curve (with a > 2) at late 
stages. However, beaming is not the only factor that leads to 
such a steep light curve. If the GRB progenitor is a massive 
star, a stellar wind environment will be created, where the 
density scales as n oc R~ 2 . Then the afterglow light curve 
will also be quite steep (Chevalier & Li 1999; Halpern et al. 
1999; Frail et al. 1999; Dai & Lu 1998a). This makes the 
problem even more complicated. 

In practical afterglow observations, we may encounter 
four kinds of optical light curves: (a) a single flat line, with 
q ~ 1.1; (b) a single steep line, with a > 2; (c) a slightly 
broken curve, with a evolves from ~ 1.1 to ~ 1.5; (d) a 
steeply broken curve, with a evolves from ~ 1.1 to > 2.1. 
let us discuss their meanings one by one: 

(a) A single flat line, with a ~ 1.1. In this case, we 
can safely say that the ejecta is not highly beamed. GRB 
970508, 971214, 980329 and 980703 may be good examples. 

(b) A single steep line, with a > 2. This may be due 
to either a highly collimated jet or an isotropic fireball in 
a wind environment. We can not distinguish them solely by 
optical light curve features. Frail et al. (1999) suggested that 
they can be distinguished by radio afterglows. GRB 980326, 
980519 and probably GRB 970228 belong to this case. 

(c) A slightly broken curve, with a evolves from 
~ 1.1 to ~ 1.5. In this case the ejecta should be isotropic, 
and the break is likely due to the relativistic-Newtonian 
transition of the remnant. 

(d) A steeply broken curve, with a evolves from 
~ 1.1 to > 2.1. In this case, the ejecta is probably highly 
collimated. The break is due to the relativistic-Newtonian 
transition of the jet. 

Showing obvious breaks in the optical light curves, GRB 
990123 and 990510 are the most hopeful candidates for 
beaming. However their late time timing index are still not 
accurately determined, we cannot definitely say that they 
belong to case (c) or (d). To determine the late time index 
accurately, we need a well defined optical light curve with 
t > 30 d. This is a difficult task. 

^From above discussions, we see that till now we have 
observed at least three kinds of optical light curves. They 
belong to case (a), case (b) and case (c) or (d) respectively. 
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